Multiplexing PCR allows the identification of within-species genetic diversity in ancient eDNA.

Sedimentary ancient DNA (sedaDNA) has rarely been used to obtain population-level data due to either a lack of taxonomic resolution for the molecular method used, limitations in the reference material or inefficient methods. Here, we present the potential of multiplexing different PCR primers to retrieve population-level genetic data from sedaDNA samples. Vaccinium uliginosum (Ericaceae) is a widespread species with a circumpolar distribution and three lineages in present-day populations. We searched 18 plastid genomes for intraspecific variable regions and developed 61 primer sets to target these. Initial multiplex PCR testing resulted in a final set of 38 primer sets. These primer sets were used to analyse 20 lake sedaDNA samples (11,200 cal. yr BP to present) from five different localities in northern Norway, the Alps and the Polar Urals. All known V. uliginosum lineages in these regions and all primer sets could be recovered from the sedaDNA data. For each sample on average 28.1 primer sets, representing 34.15 sequence variants, were recovered. All sediment samples were dominated by a single lineage, except three Alpine samples which had co-occurrence of two different lineages. Furthermore, lineage turnover was observed in the Alps and northern Norway, suggesting that present-day phylogeographical studies may overlook past genetic patterns. Multiplexing primer is a promising tool for generating population-level genetic information from sedaDNA. The relatively simple method, combined with high sensitivity, provides a scalable method which will allow researchers to track populations through time and space using environmental DNA.

The evolutionary fate of the chloroplast and nuclear rps16 genes as revealed through the sequencing and comparative analyses of four novel legume chloroplast genomes from Lupinus.

The Fabaceae family is considered as a model system for understanding chloroplast genome evolution due to the presence of extensive structural rearrangements, gene losses and localized hypermutable regions. Here, we provide sequences of four chloroplast genomes from the Lupinus genus, belonging to the underinvestigated Genistoid clade. Notably, we found in Lupinus species the functional loss of the essential rps16 gene, which was most likely replaced by the nuclear rps16 gene that encodes chloroplast and mitochondrion targeted RPS16 proteins. To study the evolutionary fate of the rps16 gene, we explored all available plant chloroplast, mitochondrial and nuclear genomes. Whereas no plant mitochondrial genomes carry an rps16 gene, many plants still have a functional nuclear and chloroplast rps16 gene. Ka/Ks ratios revealed that both chloroplast and nuclear rps16 copies were under purifying selection. However, due to the dual targeting of the nuclear rps16 gene product and the absence of a mitochondrial copy, the chloroplast gene may be lost. We also performed comparative analyses of lupine plastomes (SNPs, indels and repeat elements), identified the most variable regions and examined their phylogenetic utility. The markers identified here will help to reveal the evolutionary history of lupines, Genistoids and closely related clades.

DNA metabarcoding diet analysis for species with parapatric vs sympatric distribution: a case study on subterranean rodents.

Closely related sympatric species commonly develop different ecological strategies to avoid competition. Ctenomys minutus and C. flamarioni are subterranean rodents parapatrically distributed in the southern Brazilian coastal plain, showing a narrow sympatric zone. To gain understanding on food preferences and possible competition for food resources, we evaluated their diet composition performing DNA metabarcoding analyzes of 67 C. minutus and 100 C. flamarioni scat samples, collected along the species geographical ranges. Thirteen plant families, mainly represented by Poaceae, Araliaceae, Asteraceae and Fabaceae, were identified in the diet of C. minutus. For C. flamarioni, 10 families were recovered, with a predominance of Poaceae, Araliaceae and Asteraceae. A significant correlation between diet composition and geographical distance was detected in C. minutus, whereas the diet of C. flamarioni was quite homogeneous throughout its geographical distribution. No significant differences were observed between males and females of each species. However, differences in diet composition between species were evident according to multivariate analysis. Our results suggest some level of diet partitioning between C. flamarioni and C. minutus in the sympatric region. While the first species is more specialized on few plant items, the second showed a more varied and heterogeneous diet pattern among individuals. These differences might have been developed to avoid competition in the region of co-occurrence. Resource availability in the environment also seems to influence food choices. Our data indicate that C. minutus and C. flamarioni are generalist species, but that some preference for Poaceae, Asteraceae and Araliaceae families can be suggested for both rodents.

DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet.

Ecological understanding of the role of consumer-resource interactions in natural food webs is limited by the difficulty of accurately and efficiently determining the complex variety of food types animals have eaten in the field. We developed a method based on DNA metabarcoding multiplexing and next-generation sequencing to uncover different taxonomic groups of organisms from complex diet samples. We validated this approach on 91 faeces of a large omnivorous mammal, the brown bear, using DNA metabarcoding markers targeting the plant, vertebrate and invertebrate components of the diet. We included internal controls in the experiments and performed PCR replication for accuracy validation in postsequencing data analysis. Using our multiplexing strategy, we significantly simplified the experimental procedure and accurately and concurrently identified different prey DNA corresponding to the targeted taxonomic groups, with ≥ 60% of taxa of all diet components identified to genus/species level. The systematic application of internal controls and replication was a useful and simple way to evaluate the performance of our experimental procedure, standardize the selection of sequence filtering parameters for each marker data and validate the accuracy of the results. Our general approach can be adapted to the analysis of dietary samples of various predator species in different ecosystems, for a number of conservation and ecological applications entailing large-scale population level diet assessment through cost-effective screening of multiple DNA metabarcodes, and the detection of fine dietary variation among samples or individuals and of rare food items.

DNA from soil mirrors plant taxonomic and growth form diversity.

Ecosystems across the globe are threatened by climate change and human activities. New rapid survey approaches for monitoring biodiversity would greatly advance assessment and understanding of these threats. Taking advantage of next-generation DNA sequencing, we tested an approach we call metabarcoding: high-throughput and simultaneous taxa identification based on a very short (usually <100 base pairs) but informative DNA fragment. Short DNA fragments allow the use of degraded DNA from environmental samples. All analyses included amplification using plant-specific versatile primers, sequencing and estimation of taxonomic diversity. We tested in three steps whether degraded DNA from dead material in soil has the potential of efficiently assessing biodiversity in different biomes. First, soil DNA from eight boreal plant communities located in two different vegetation types (meadow and heath) was amplified. Plant diversity detected from boreal soil was highly consistent with plant taxonomic and growth form diversity estimated from conventional above-ground surveys. Second, we assessed DNA persistence using samples from formerly cultivated soils in temperate environments. We found that the number of crop DNA sequences retrieved strongly varied with years since last cultivation, and crop sequences were absent from nearby, uncultivated plots. Third, we assessed the universal applicability of DNA metabarcoding using soil samples from tropical environments: a large proportion of species and families from the study site were efficiently recovered. The results open unprecedented opportunities for large-scale DNA-based biodiversity studies across a range of taxonomic groups using standardized metabarcoding approaches.

Using next-generation sequencing for molecular reconstruction of past Arctic vegetation and climate.

Palaeoenvironments and former climates are typically inferred from pollen and macrofossil records. This approach is time-consuming and suffers from low taxonomic resolution and biased taxon sampling. Here, we test an alternative DNA-based approach utilizing the P6 loop in the chloroplast trnL (UAA) intron; a short (13-158 bp) and variable region with highly conserved flanking sequences. For taxonomic reference, a whole trnL intron sequence database was constructed from recently collected material of 842 species, representing all widespread and/or ecologically important taxa of the species-poor arctic flora. The P6 loop alone allowed identification of all families, most genera (>75%) and one-third of the species, thus providing much higher taxonomic resolution than pollen records. The suitability of the P6 loop for analysis of samples containing degraded ancient DNA from a mixture of species is demonstrated by high-throughput parallel pyrosequencing of permafrost-preserved DNA and reconstruction of two plant communities from the last glacial period. Our approach opens new possibilities for DNA-based assessment of ancient as well as modern biodiversity of many groups of organisms using environmental samples.

Assessment of microbial communities by graph partitioning in a study of soil fungi in two Alpine meadows.

Understanding how microbial community structure and diversity respond to environmental conditions is one of the main challenges in environmental microbiology. However, there is often confusion between determining the phylogenetic structure of microbial communities and assessing the distribution and diversity of molecular operational taxonomic units (MOTUs) in these communities. This has led to the use of sequence analysis tools such as multiple alignments and hierarchical clustering that are not adapted to the analysis of large and diverse data sets and not always justified for characterization of MOTUs. Here, we developed an approach combining a pairwise alignment algorithm and graph partitioning by using MCL (Markov clustering) in order to generate discrete groups for nuclear large-subunit rRNA gene and internal transcript spacer 1 sequence data sets obtained from a yearly monitoring study of two spatially close but ecologically contrasting alpine soils (namely, early and late snowmelt locations). We compared MCL with a classical single-linkage method (Ccomps) and showed that MCL reduced bias such as the chaining effect. Using MCL, we characterized fungal communities in early and late snowmelt locations. We found contrasting distributions of MOTUs in the two soils, suggesting that there is a high level of habitat filtering in the assembly of alpine soil fungal communities. However, few MOTUs were specific to one location.

Origin and fate of repeats in bacteria.

We investigated 53 complete bacterial chromosomes for intrachromosomal repeats. In previous studies on eukaryote chromosomes, we proposed a model for the dynamics of repeats based on the continuous genesis of tandem repeats, followed by an active process of high deletion rate, counteracted by rearrangement events that may prevent the repeats from being deleted. The present study of long repeats in the genomes of Bacteria and Archaea suggests that our model of interspersed repeats dynamics may apply to them. Thus the duplication process might be a consequence of very ancient mechanisms shared by all three domains. Moreover, we show that there is a strong negative correlation between nucleotide composition bias and the repeat density of genomes. We hypothesise that in highly biased genomes, non-duplicated small repeats arise more frequently by random effects and are used as primers for duplication mechanisms, leading to a higher density of large repeats.

Study of intrachromosomal duplications among the eukaryote genomes.

Complete eukaryote chromosomes were investigated for intrachromosomal duplications of nucleotide sequences. The analysis was performed by looking for nonexact repeats on two complete genomes, Saccharomyces cerevisiae and Caenorhabditis elegans, and four partial ones, Drosophila melanogaster, Plasmodium falciparum, Arabidopsis thaliana, and Homo sapiens. Through this analysis, we show that all eukaryote chromosomes exhibit similar characteristics for their intrachromosomal repeats, suggesting similar dynamics: many direct repeats have their two copies physically close together, and these close direct repeats are more similar and shorter than the other repeats. On the contrary, there are almost no close inverted repeats. These results support a model for the dynamics of duplication. This model is based on a continuous genesis of tandem repeats and implies that most of the distant and inverted repeats originate from these tandem repeats by further chromosomal rearrangements (insertions, inversions, and deletions). Remnants of these predicted rearrangements have been brought out through fine analysis of the chromosome sequence. Despite these dynamics, shared by all eukaryotes, each genome exhibits its own style of intrachromosomal duplication: the density of repeated elements is similar in all chromosomes issued from the same genome, but is different between species. This density was further related to the relative rates of duplication, deletion, and mutation proper to each species. One should notice that the density of repeats in the X chromosome of C. elegans is much lower than in the autosomes of that organism, suggesting that the exchange between homologous chromosomes is important in the duplication process.

Analysis of intrachromosomal duplications in yeast Saccharomyces cerevisiae: a possible model for their origin.

The complete genome of the yeast Saccharomyces cerevisiae was investigated for intrachromosomal duplications at the level of nucleotide sequences. The analysis was performed by looking for long approximate repeats (from 30 to 3,885 bp) present on each of the chromosomes. We show that direct and inverted repeats exhibit very different characteristics: the two copies of direct repeats are more similar and longer than those of inverted repeats. Furthermore, contrary to the inverted repeats, a large majority of direct repeats appear to be closely spaced. The distance (delta) between the two copies is generally smaller than 1 kb. Further analysis of these "close direct repeats" shows a negative correlation between delta and the percentage of identity between the two copies, and a positive correlation between delta and repeat length. Moreover, contrary to the other categories of repeats, close direct repeats are mostly located within coding sequences (CDSs). We propose two hypotheses in order to interpret these observations: first, the deletion/conversion rate is negatively correlated with delta; second, there exists an active duplication mechanism which continuously creates close direct repeats, the other intrachromosomal repeats being the result, by chromosomal rearrangements of these "primary repeats."

A comparative study of duplications in bacteria and eukaryotes: the importance of telomeres.

The genomes of three bacteria (Haemophilus influenzae, Mycoplasma genitalium, and Escherichia coli) and two eukaryotes (Saccharomyces cerevisiae and Caenorhabditis elegans) were compared. The distribution of their putative open reading frames (ORFs) was studied, and several conclusions were drawn: (1) All of these genomes, even the smallest, exhibit a significant proportion (7%-30%) of duplicated ORFs. This proportion is a function of genome size and appears unrelated to the bacteria/eukaryote division. (2) Some of these ORFs constitute families of up 20 or more members. (3) The levels of sequence similarity within these families are highly variable and their distribution is different among bacteria and eukaryotes. (4) In yeast, there are topological relationships between members of the same family. The paired ORFs are frequently in the same orientation with regard to their respective telomeres and located at comparable distances from them.

Unsuspected prophage-like elements in Salmonella typhimurium.

We present evidence for the existence of two large (approximately 50 kb) excisable segments in the chromosome of Salmonella typhimurium. The two elements--designated Gifsy-1 and Gifsy-2--cover, respectively, the 57 units and the 24 units of the genetic map where they contribute indicative rare restriction sites. The two elements are closely interrelated and both contain a region of sequence similarity to the recE locus of the Rac prophage of Escherichia coli. Mutations within this region of Gifsy-1 yield the classical 'Sbc' phenotype: they suppress the recombination defect of recB mutants, apparently by activating a normally silent recE-like gene. At the same time, these 'sbcE' mutations activate a Xis-type function that promotes excision of one or other of the two elements. Predictably, curing of Gifsy-1 results in the loss of recB mutant suppression. Surprisingly, the suppressor phenotype is also lost in cells cured for Gifsy-2 even though the Gifsy-1-associated sbcE mutation is still present. Moreover, the excision frequency of Gifsy-1 drops dramatically in Gifsy-2-cured cells. Thus, both elements must co-operate in the activation of recombination and excision functions. Overall, the data presented here suggest that Gifsy-1 and Gifsy-2 are cryptic prophages. They are distinct from previously described Fels prophages. Unlike Fels, they are not specific to S. typhimurium strain LT2 since they are both also found in a virulent S. typhimurium isolate (ATCC 14028s).

The nucleotide sequence of Saccharomyces cerevisiae chromosome VII.

The complete nucleotide sequence of Saccharomyces cerevisiae chromosome VII has 572 predicted open reading frames (ORFs), of which 341 are new. No correlation was found between G+C content and gene density along the chromosome, and their variations are random. Of the ORFs, 17% show high similarity to human proteins. Almost half of the ORFs could be classified in functional categories, and there is a slight increase in the number of transcription (7.0%) and translation (5.2%) factors when compared with the complete S. cerevisiae genome. Accurate verification procedures demonstrate that there are less than two errors per 10,000 base pairs in the published sequence.

Sequence of a 39,411 bp DNA fragment covering the left end of chromosome VII of Saccharomyces cerevisiae.

We have sequenced a DNA fragment of 39,411 bp which includes part of the left telomere of chromosome VII of Saccharomyces cerevisiae. We have identified 19 open reading frames (ORFs); six correspond to known yeast genes (ADH4, FZF1, HKB, RTG2, HFM1 and PDE1), nine have similarity with other genes and four exhibit no significant similarity with any known gene. The average size of these ORFs seems to be related to their location, the eight ORF's nearest the telomere being shorter than the 11 others. These two groups of genes are separated by a region of 4.5 kb devoid of significant ORFs. One ORF, NRF120, is a new member of the seripauperine family, represented once in all sequenced yeast chromosomes, in a subtelomeric location.

Cloning of a large gene cluster involved in Erwinia amylovora CFBP1430 virulence.

Phage MudIIPR13 insertional mutagenesis of Erwinia amylovora CFBP1430 allowed us to isolate 6900 independent CmR mutants. The frequencies of different auxotrophs in this population indicated that MudIIPR13 had inserted randomly in E. amylovora. Screening of 3500 CmR mutants on (i) apple calli and (ii) pear and apple seedlings led to the isolation of 19 non-pathogenic prototrophic single mutants, four of which expressed a LacZ+ hybrid protein. Expression of the fusion proteins was temperature sensitive. The 19 mutants could be separated into two classes according to their behaviour on tobacco: 13 were unable to elicit the hypersensitive response on tobacco (Hrp-) while six still could (Dsp-). The 19 MudIIPR13 insertions all mapped in the same virulence region. The MudIIPR13 insertions of Hrp- mutants were all clustered on the left part of this region, while the MudIIPR13 insertions of Dsp- mutants were located on the right part. All of the mutants except one, which proved to have a large deletion of the entire virulence region, could be complemented functionally by cosmids from an E. amylovora CFBP1430 genomic library. No hybridization was observed between the cosmid pPV130, which complemented 12 hrp::MudIIPR13 mutations, and the hrp genes from Pseudomonas syringae pv. phaseolicola (Lindgren et al., 1986), P. syringae pv. tomato (N.J. Panopoulos, unpublished data) or P. solanacearum (Boucher et al., 1987). Further analysis of the large virulence region will allow mapping of the border of the virulence region and facilitate the study of the function and regulation of the hrp and dsp genes.